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ABSTRACT 



u 

^' 

Aims. We explore in this paper the ability of spatially resolved spectroscopic measurements of the SZ effect (SZE) to deter- 
^ , mine the temperature profile of galaxy clusters. We derive a general formalism for the thermal SZE in galaxy clusters with a 
• non-uniform temperature profile that can be applied to both cool-core clusters and non-cool core cluster with an isothermal or 
C/3 I non-isothermal temperature structure. 

. Methods. We derive an inversion technique through which the electron distribution function can be extracted from spectroscopic 
■ SZE observations over a wide frequency range. We study the fitting procedure to extract the cluster temperature from a set of 
I simulated spatially resolved spectroscopic SZE observations in diff'erent bands of the spectrum, from 100 to 450 GHz. 
^ Results. The results of our analysis for three different cluster prototypes (A2199 with a low-temperature cool core, Perseus with 
a relatively high-temperature cool core, Ophiuchus with an isothermal temperature distribution) provide both the required 
^SJ , precision of the SZE observations and the optimal frequency bands for a determination of the cluster temperature similar or 
CNl ■ better than that obtainable from X-ray observations. The precision of SZE-derived temperature is also discussed for the outer 
regions of clusters. We also study the possibility to extract, from our method, the parameters characterizing the non-thermal 
. SZE spectrum of the relativistic plasma contained in the lobes of radio galaxies as well as the spectrum of relativistic electrons 
y—{ ' co-spatially distributed with the thermal plasma in clusters with non-thermal phenomena. 

^\ Conclusions. We find that the next generation SZE experiments with spectroscopic capabilities, like those using FTS spectrome- 
. ters with imaging capabilities, can provide precise temperature distribution measurements over a large range of radial distances 
' for galaxy clusters even out to substantial redshifts. 
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1. Introduction ies provide information on the thermodynamic state of 

• • „ 1 , the ICM, allow to measure the total mass of these sys- 

Galaxy clusters are the largest eravitationally bound con- , i i j j. j_- -i-i, ■ i- /- 

, , , tems through the hydrostatic equilibrium equation (see 

tamer ol dinuse baryons m the universe. These baryons ac- . onnc; ri i f p j r,oo'7^ j 

, , , 1 e.g. Arnaud 2005, Colatrancesco & Giordano 2007), and 

comodate m the cluster gravitational potential well domi- n.. ^ ■ ^ 4.1, i ^ i-i 

ii^iiT^inT 11 - i- i 1 allow to set constraints on the relevant cosmological pa- 

nated by Cold Dark Matter and show, m their hotter phase . /r-i r-. \ ^^ 1 ii _l 1 r ii 1 _l 

' Tr^^^r -^u rameters (S2m,S2D£;, Sib) through the study of the cluster 

(i.e. the hot Intra Cluster Medium - hcrcaltcr ILM - with , ,. , ,. 

^ 7 s N 1 population evolution, 
temperature ~ 10 — 10 K), a complex temperature 

distribution in the cluster atmospheres as indicated by ob- On the cosmological side, once the gas mass is deter- 
servations (see, e.g., Arnaud 2005 for a review; Pratt et mined from X-ray observations in the deprojection or fit- 
al. 2007 for recent XMM-Newton results, Ehlert & Ulmer ting analysis (see, e.g., Fabian et al. 1981; Ettori & Fabian 
2009 for recent Chandra results, Sato et al. 2009 for recent 1999) and the total mass is estimated through the hydro- 
Suzaku results) and numerical simulations (see Borgani et static equilibrium condition, it is possible to derive the 
al. 2008 for a recent review). cluster gas mass fraction that can be used to constrain 
Spatially resolved studies of the ICM temperature dis- the cosmological density parameter, if combined with pri- 
tribution in galaxy clusters are crucial for both astro- mordial nucleo-synthesis calculations (e.g.. White et al. 
physical and cosmological apphcations since these stud- 1993, Ettori et al. 2009). Cosmological apphcations of 
the cluster population evolution require precise measure- 
Send offprint requests to: S. Colatrancesco ments of the cluster temperature and of the correlation be- 
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tween temperature and other relevant physical quantities, 
in order to describe correctly the cluster population (to- 
tal mass, luminosity, Compton parameter, number counts, 
redshift distribution, luminosity function). If well cali- 
brated, the slope and evolution of cluster scaling relations, 
such as gas mass versus temperature (e.g., Colafrancesco 
& Vittorio 1994, Colafrancesco et al. 1994, 1997, Voit 
2000, Finoguenov, Reiprich, & Bochringer 2001) and clus- 
ter size versus temperature (e.g., Colafrancesco, et al. 
1997, Mohr et al. 2000; Verde et al. 2001), can also be 
used to constrain cosmological and structure formation 
models. 

On the astrophysical side, the precise determination 
of important quantities such as entropy, pressure, and - 
under the assumption of hydrostatic equilibrium - the to- 
tal cluster mass, is dependent on accurate estimation of 
the temperature profile. Accurate temperature profiles are 
fundamental in determining the gas entropy distribution 
(see Lloyd-Davics, Ponman, & Cannon 2000 and refer- 
ences therein) , which is a powerful tool to explore the role 
of non-gravitational processes that could alter the specific 
thermal energy in the ICM (Colafrancesco 2008a for a re- 
view, Kang et al. 2007, Colafrancesco & Giordano 2006, 
2007), to assess the structure of cool cores and their sur- 
rounding atmosphere in order to probe the interplay of 
non-thermal and thermal particle distributions in clusters 
(e.g., Colafrancesco, Dar & DeRujula 2004, Guo & Ho 
2008, Colafrancesco & Marchegiani 2008 and references 
therein) and to assess the amount of non-thermal (point- 
like and/or diffuse) emission mechanisms detectable in 
the hard X-ray domain (see discussion in Colafrancesco 
& Marchegiani 2009). 

The standard methodology to recover the ICM temper- 
ature distribution makes use of high resolution spectro- 
scopic X-ray observations (XMM, Chandra, Suzaku and 
previously ROSAT, ASCA, and BeppoSAX). Such obser- 
vations provide the key measurable characteristics of the 
ICM, i.e., the temperature and density of this plasma. 
Because of limited photon statistics it is usual to measure 
the density and temperature in terms of radial profiles 
(see, e.g., Pratt et al. 2007). However, while the density 
of the ICM is relatively easy to measure from the surface 
brightness profile of a galaxy cluster, precise temperature 
determination requires high photon statistics to build, and 
fit, a spectrum. Therefore. ICM temperature profiles are 
typically determined with considerably less spatial reso- 
lution than density profiles. The measurement of radial 
temperature profiles is further complicated by tlic^ den- 
sity squared dependence of the thermal bremsstrahlung 
X-ray emissivity, Sbrem oc nl{r)Te^'^{r). The steep drop 
of the X-ray surface brightness with distance from the 
centre, combined with the background from cosmic, solar 
and instrumental sources, makes, hence, accurate X-ray 
measurement of the temperature distribution at large dis- 
tances from the cluster centre a technically challenging 
task (see discussion in Pratt et al. 2007). 



ICM temperature can also be measured, alternatively, 
by using the Inverse Compton sc;attc;ring of CMB pho- 
tons off thermal electrons residing in the cluster atmo- 
sphere - the Sunyaev-Zel'dovich effect (SZE; see Sunyaev 
& Zel'dovich 1980, Birkinshaw 1999, Colafrancesco 2007 
for reviews). This effect provides a CMB temperature 
change 

AT ctt f 

= / Mn.kBT, ■ g{x) (1) 

To meC^ Ji 

that depends on the cluster temperature directly from 
its amplitude rx y = MneksTe (where ax is the 

Thomson cross section, nie is the electron mass, Ug is the 
electron density, Tg is the cluster temperature, the integral 
is performed along the line of sight £), and from its spec- 
trum oc g{x) (where x = hu/ksTo is the frequency normal- 
ized to the CMB energy, h is the Planck constant, fcs is the 
Boltzmann constant and Tq is the present-day CMB tem- 
perature). The temperature dependence is explicitly in- 
cluded in the function g{x) through the relativistic effects 
of the Compton Scattering that are more prominent in 
the high frequency region of the SZE spectrum at z/ ^ 300 
GHz (sec Birkinshaw 1999, Colafrancesco, Marchegiani & 
Palladino 2003, Colafrancesco 2007 for a review), where it 
is therefore possible to measure directly the ICM temper- 
ature. This is evident looking at the different temperature 
dependence of the three basic spectral features that char- 
acterize the thermal SZE signal: 

i) the minimum in its intensity located at the frequency 

xth,min « 2.265(1 - 0.0927f?e + 2.386*^) 

-Ft(-0.00674 + 0.4666'e) , (2) 

where 9e = ksTe/meC^, whose value depends weakly on 
the electron spectrum (i.e. on Tg and rig) and equals ~ 
2.26; 

ii) the crossover frequency, xq, whose value depends on 
the electron pressure / energy density and optical depth 

xthfi ~ a{Te) + rb{Te) (3) 

with a(Te) = 3.830(1 + 1.1620^ - 0.81446(2) and b{Te) = 
3.0216»e - 8.6726'2, and is found at a frequency > 3.83 for 
increasing values of Tg (the value xth.o = 3.83 is found in 
the non-relativistic limit, or in the limit Tg ^0); 

iii) the maximum of its intensity located at the frequency 

Xth,max « 6.511(1 + 2.410e-4.960g) 

-FT(0.0161-|-8.166»e- 35.96»^) (4) 

and depends sensitively on the nature of the elec- 
tron population and on its energy (momentum) spec- 
trum (see Dolgov et al. 2001, Colafrancesco et al. 2003, 
Colafrancesco et al. 2009 for the case of electrons with a 
thermal spectrum; see also Colafrancesco et al. 2003 and 
Colafrancesco 2004, Colafrancesco 2005, Colafrancesco 
2007, Colafrancesco 2008b for the case of electron with 
different spectra). 

The SZE can be used, therefore, as an alternative 
probe of the ICM temperature provided that detailed 
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spectral measurements extending out to high frequencies 
can be obtained (see Colafrancesco 2007 for a review, see 
Colafrancesco et al. 2009). 

The linear dependence of the SZE temperature change AT 
from the electronic density Ue and temperature (see eq. 
[1]) makes it possible - in addition - to enhance the sen- 
sitivity of the SZE measurements in the outer regions of 
the cluster while maintaining an adequate sensitivity also 
in the inner parts. Therefore SZE observations could pro- 
vide a more uniform and spatially extended coverage of 
the temperature profile in galaxy clusters than X-ray ob- 
servations. 

Little or nothing is known about cluster tempera- 
ture profiles obtained directly from SZE observations. 
Present-day and planned SZE observations from both 
ground based experiments (multi-band bolometers and in- 
terferometers) or coming space-borne experiments (e.g., 
OLIMPO and PLANCK) require complementary X-ray 
observations to put tight constraints on the cluster tem- 
perature, even though a somewhat more moderate pre- 
cision on cluster temperature seems to be sufficient to 
achieve a reasonable statistical knowledge on the CDM 
power spectrum, namely on erg (Juin et al. 2007). 
We must also stress, in this context, that the use of X-ray 
observed quantities in the SZE data analysis is, however, a 
delicate procedure, because the SZE-derived temperature 
is weighted by the Compton parameter y. It is known, for 
instance, that the derived peculiar velocity is systemati- 
cally shifted by ^ 10 — 20% due to this reason (see e.g., 
Hansen 2004). 

While observations of the SZE are becoming increasingly 
accurate (e.g., LaRoque et al. 2002; De Petris et al. 2002; 
Battistelli et al. 2003, Halverson et al. 2009, Staniszewski 
et al. 2009), most analyses of SZ data are still made un- 
der the simplifying assumption of iso-thermality, or single 
temperature plasma. Several studies have also considered 
the ability of future SZ observations to extract the clus- 
ter physical parameters (see, e.g., Knox, Holder & Church 
2004; Aghanim, Hansen & Lagache 2005), but always un- 
der the assumption of iso-thermality. 
A desirable solution with the SZE observations is to 
achieve simultaneous good angular resolution and high 
sensitivity in order to deproject the SZE signal properly 
throughout the whole cluster atmosphere. 

In this paper we discuss specifically the strategy that 
will be more effective for obtaining spatially-resolved tem- 
perature measurements which are only based on SZE ob- 
servations. Such strategy needs to overcome the following 
problems: i) achieve spatially resolved spectral measure- 
ments of the cluster temperature; ii) measure the temper- 
ature profile over a wide spatial range in order to disentan- 
gle the cool cores region from the outer cluster atmosphere 
(or regions with different temperature in the overall atmo- 
sphere) ; iii) obtain detailed and independent temperature 
measurements which require a wide frequency band cov- 
erage with sufficient spectral resolution. These points are 
the main issues that our study wants to address. 



The outline of the paper is the following. We describe 
the general theory of the SZE in Sect. 2, and we general- 
ize this theory to the case of clusters with non-uniform 
T-profiles, like e.g. cool core clusters, in Sect. 3. We de- 
scribe a fitting procedure for spatially resolved spectro- 
scopic SZE observations in Sect. 4, where we first dis- 
cuss the case of the thermal SZE. The analysis of the 
non-thermal SZE and the combination of thermal and 
non-thermal SZE are also discussed, for completeness, in 
Sect. 5. The non-thermal SZE is important both for the 
intrinsic study of non-thermal and/or relativistic plasmas 
in clusters and for determining their role as biases for the 
study of the thermal plasma (i.e. the ICM). We present in 
Sect. 6 the application of our study for two specific galaxy 
clusters: Perseus, which has a cool core, and Ophiuchus, 
which is approximately isothermal. We discuss our results 
and summarize our conclusions in the final Sect. 7. 

Throughout the paper, we use a flat, vacuum- 
dominated cosmological model with ^lm = 0.3, f^A = 0.7 
and Hq ^ 70 km s^^ Mpc^^. 

2. The SZE: theory 

In order to derive the expression of the SZE which is valid 
for a galaxy cluster with spatially varying density and tem- 
perature, we need to recall the general expressions for the 
SZE formulae. 

We use here the formalism presented by Wright (1979), 
and subsequently developed by Birkinshaw (1999), Enfilin 
& Kaiser (2000), Colafrancesco, Marchegiani & Palladino 
(2003). We note that Bcehm & Lavalle (2009) recently 
proposed a covariant formalism for the SZE, claiming, er- 
roneously, that the formalism presented in Wright (1979) 
and further developed by other authors was incorrect. We 
have verified (both numerically and analytically) that the 
two approaches, in the Thomson limit, are fully equivalent 
and provide the same results (see Appendix A for details). 
In addition to our verification, Nozawa & Kohyama (2009) 
have also shown analytically that the two approaches are 
fully equivalent, pointing out the errors made by Bcehm 
& Lavalle (2009). 

The analytic expression for the upscattered CMB spec- 
trum is: 

/ + 00 
/o(xe-")P(s)ds (5) 

(see e.g. Colafrancesco, Marchegiani & Palladino 2003) 
where 

with X = hv/kBT{). 

The function P{s) is given by the following expression 

n=0 
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where 
T = ax I UeM 



is the electron optical depth and 



(8) 



(9) 



n times 



where the symbol indicates the convolution product, 
and 



fe{p)Ps{s,p)dp 



(10) 



The function fe{p) is the electron momentum distribu- 
tion function (with p = f3j) and is normalized as to give 
/o°° fe{p)dp = 1. The function Ps{s,p) is known from the 
basic physics of the Compton scattering and is given by 



Ps{s,p)ds = P{e'',p)e'ds 
with 

Z\l-t 



(11) 



P{t,p) 



+ 



i2pH 
3(1 + <) 



1 + (10 + + 4p'')i + t'\ + 

3 + 3p2+p4 



3_+2p2 
2p 



(2arcsinh(p) - |ln(t)|) 



(12) 



Our aim is to invert eq.(l5]) in order to derive the spectral 
and spatial parameters of the electrons distribution from 
measurements of the SZE only. We derive here an analyt- 
ical solution to this problem, while in Sect. 4 we present 
a numerical procedure to solve it. 

The approximated expression of A/ (a;) at first order 
in T is 



where 



lQ{xe-'')Pi{s)ds . 



Thus, by using eq. fTU)) . we can write 
A/(a;) 



+ Io{x)= / k{x,p)feip)dp 
T Jo 

where 

/+00 
Io{xe~'')Ps{s,p)ds . 
-oo 

Introducing the following definitions 



nux) 



k{x,p)fe{p)dp 



I(x) 

Ir{x) ^ ±l+Io{x), 



(13) 
(14) 

(15) 

(16) 

(17) 
(18) 



the previous eg 1 151 writes as 

ir{x)^nfe{x) . (19) 

Its formal solution is given by 

f,{x)^n-^ir{x) (20) 

with 

n-^H = I, (21) 
where T is the identity operator. 

3. The SZE for a galaxy cluster with a radial 
temperature profile 

In order to apply the general formalism described in the 
previous Section to a cluster with a spatially varying tem- 
perature and density profiles (i.e., a cluster with a cool 
core or a non-isothermal cluster), we must consider that: 

— equations ©-([1]) still hold. In particular, the expres- 
sion for the optical depth remains the same because it 
depends only on the gas density profile and not on the 
cluster temperature; 

— the expression for the function -P(s), and its series ex- 
pansion in terms of powers of r, remains formally un- 
changed, because it gives the probability that a CMB 
photon suffers a (logarithmic) change in the frequency 
s = ln{v' /v) by traversing the whole cluster; 

— the expression for the function Pi{s) (i.e. eg. fTU]) must 
be changed. In fact, the expression given in eq.([TO]) as- 
sumes that the electron momentum distribution func- 
tion /e(p) is the same throughout the cluster atmo- 
sphere. Such an assumption is not correct for a cluster 
with a cool core (or a non-isothermal cluster) since its 
temperature, and hence the electron momentum dis- 
tribution, changes with the cluster radius. 

We must, therefore, find an appropriate way to calcu- 
late the function Pi{s) taking into account the fact that 
the CMB photon, while traversing the cluster atmosphere, 
finds regions with different temperature and density. In 
order to calculate the redistribution function Pi{s) ap- 
propriately, we can calculate - along the line of sight £ 
- the local electron momentum distributions weighted by 
its local density. 

Describing the thermal electron density as 



(22) 



it is possible to derive the average density-weighted mo- 
mentum distribution along the line of sight ^ as 



Up) 



where 



JJeip;T,{r))g,{r)de 
ji9e{r)dl 



(23) 



(24) 
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with r] = (meC^)/(/esT'e) a-nd K2{r]) being the modified 
Bessel function of second kind (see, e.g., Abramowitz & 
Stegun 1965), which ensures the correct normahzation 



/•oo 

Jo 



(25) 



The function /e also satisfies the normalization given in 
eg. ([^5]) . as can be verified by the foUowing equation 



fe{p)dp = 



/o°° d.pf^fe{p;Te{r))ge{r)dl 

j^ge{r)dl 
Ji 9e {r)d£ dpfe (p; Te{r)) 



Ji9e{r)de 



1 



(26) 



where we use eq. (^5]) . 

We can, therefore, re- write ea. pH)) for a cluster with 
a generic temperature distribution by using the averaged 
electron momentum distribution 



fe{p)Ps{s,p)dp , 



(27) 



where fe{p) is given in eq. (l23|) . This function allows us 
to calculate the frequency redistribution probability for 
multiple scattering and, hence, the analytic expression of 
the SZE up to the required approximation order by using 
eqs. d!])-®. 

4. Fitting procedure: the thermal SZE 



Eq. (l20|) cannot be easily solved by analytically methods. 
Therefore, we describe in this section a numerical proce- 
dure by which it is possible to derive the cluster temper- 
ature profile starting from observations of the SZE with 
appropriate spatial and spectral resolution. 
The quantitative analysis that we present here refers to 
the specific case of the cluster A2199 which has a cool 
core. 

The procedure that we follow consists of simulating 
the thermal SZE observation of the cluster A2199 with pa- 
rameters (density and temperature) taken from the avail- 
able X-ray observations, and then extract the best fit val- 
ues of the ICM temperature and optical depth from the 
simulated SZE observations. Note that in this section we 
work, for convenience, with the first-order approximation 
to the thermal SZE, and in Sect. 5 we will discuss the con- 
sequences of using this approximation. The temperature 
profile obtained from SZE measurements taken at vari- 
ous cluster radii is then de-projected under the constraint 
(prior) that the density profile derived from SZE observa- 
tions is consistent with the density profile obtained from 
X-ray measurements: we use this procedure (which hence 
assumes the X-ray density prior) in order to have a fast 
enough computing time that allows us to explore the de- 
tails of our study. Note, however, that the density con- 
straint is not necessary, in principle, because the cluster 
density profile can be obtained, self-consistently, from the 



optical depth profile measured from SZE measurements. 
This fully consistent procedure requires, however, much 
longer computing times since the number of free parame- 
ters to be fitted to the observations is larger. 

For the reference case of A2199 that we want to dis- 
cuss in this section, our analysis goes thorugh the following 
steps: i) we have taken X-ray information on A2199 from 
Chandra observation (Johnstone et al. 2002); ii) using the 
density and temperature profiles derived from X-ray ob- 
servations we calculated (using the procedure described 
in Sect. 3 above) the SZE spectrum at eight different pro- 



jected radii Vp from the cluster center, i.e. 



0, 5, 10, 



15, 20, 30, 50 and 100 kpc; iii) for each radial bin, we have 
sampled the SZE spectrum at six different frequencies, i.e. 
V = 300, 320, 340, 360, 380 and 400 GHz (see Fig. [I]); iv) 




Fig. 1. The SZE spectrum obtained for the cluster A2199 at 
different projected radii is shown together with the frequency 
sampling derived, for each curve, with a 0.1% uncertainty. 
From top to bottom curves refer to projected radii of 0, 5, 
10, 15, 20, 30, 50, 100 kpc. Note that spectra for 10 and 15 
kpc are almost superposed due to the fact that beyond this 
radial distance the IC gas temperature measured by Chandra 
decreases with increasing radius (Johnstone et al. 2002). 



for each radius, we fitted the six frequency SZE data with 
a relativistic model of the thermal SZE leaving free two 
parameters, i.e. fc^Te and r, and assuming for each ex- 
perimental data point an uncertainty of 0.1%. This fitting 
procedure yields, for each projected radius, a value of the 
projected temperature and optical depth (see Fig. [J); v) 
in order to deproject the fitted temperature profile, we 
have assumed a 3-D theoretical temperature profile of the 
form: 



Teij) — Tint + {Text — Tint) 



{r/rc 



1 + {rlrcY 



(28) 



6 



S. Colafrancesco and P. Marchegiani: SZ effect and cluster temperature profiles 



4 - 



20 40 60 80 100 120 
fp (kpc) 



O.OlOf 



0.008 ■ 



0.006 ■ 



0.004 - 



0.002 - 



0.000 [ 



60 80 



1 00 1 20 



r (kpc) 



Fig. 2. The radial profile of the projected temperature (top 
panel) and optical depth (bottom panel) as derived from the fit 
to the SZE spectra using data sampled in the frequency range 
300-400 GHz with uncertainty of 0.1% (see FigO. 

(see e.g. Piffaretti et al. 2005). To decrease the number 
of free parameters in the fit we have assumed the ICM 
density profile ne(r) as obtained from X-ray data. In the 
fitting procedure, the projected temperature is calculated 
by weighting the temperature profile with the gas density 
and integrating along the los: 



Te,p{Tp) 



2/«n,(r)re(r)^ 



fdr 



=dr 



(29) 



where R is the cluster radius (and it is i? = 200 kpc for 
the case of A2199). Assuming the value of Text from X- 
ray observations, the number of free parameters reduces 



to three. In this way it is possible to fit the projected 
temperature as obtained from the simulated SZE obser- 
vations and then derive the temperature profile parame- 
ter (in eqUS]); assuming ksText = 4.2 keV for the case of 
A2199, the best fit parameters are ksTint — 0.0 ±0.5 keV, 
Tc = 4.9 ± 0.6 kpc and /.t = 1.0 ± 0.1. The best-fit temper- 
ature profile curve is shown in FigOand it is compared to 
the data on the de-projected temperature as derived from 
X-ray observations. 




-20 20 40 60 80 100 120 
r (kpc) 



Fig. 3. The deprojected temperature profile curve (solid) de- 
rived from SZE spectral data is compared with the deprojected 
temperature data as obtained from X-ray measurements of 
A2199 (Johnstone et al. 2002). 

The same procedure has been also repeated assum- 
ing that the SZE observation sampled from Fig. [T] has an 
uncertainty of 1%. The radial profile of the projected tem- 
perature obtained in this case is shown in Fig. 21 In this 
case, the de-projection of the temperature profile using 
the template given in eq. (j28p does not provide, however, 
an acceptable fit. For this reason we conclude that an un- 
certainty of order of 1% in the SZE data is not sufficient to 
derive a detailed profile of the ICM temperature in clus- 
ters like A2199. 

The whole procedure has been repeated sampling the 
SZE observations in the low- frequency range 100-200 GHz 
{v = 100, 120, 140, 160, 180, 200 GHz). In such a case, we 
find that, even assuming uncertainties of 0.1%, it is not 
possible to reproduce in a satisfactory way the radial pro- 
file of the projected temperature as derived from the fit by 
using the model in eq. ((28)) (see Figl5]). This is because the 
shape of the SZE spectrum at low frequencies is less sensi- 
tive to the cluster temperature than the shape of the high- 
frequency part of the SZE spectrum. Therefore, in order 
to obtain detailed information on the cluster temperature 
profile one must use high-frequency spectral observations 
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5. Fitting procedure: the non-thermal SZE 



4 - 



20 40 



60 
(kpc) 



80 100 120 



Fig. 4. The projected temperature profile obtained from the 
fit to SZE spectra using data points sampled in the 300-400 
GHz range and assuming an uncertainty of 1%. 



in the range 300-400 GHz. This is a crucial requirement to 
plan an observational strategy of SZE observations with 
the goal of deriving strong constraints to the fundamental 
parameters of the cluster. 



2 r 



OL 



40 60 
r (kpc) 



Fig. 5. The projected temperature radial profile obtained from 
the fit to the SZE spectral data by using data points sampled 
in 100-200 GHz range and assuming uncertainties of 0.1%. 



Here we consider the effects of the variation of the param- 
eters of the non-thermal electron spectrum on the relative 
non-thermal SZE. 

First, we study the spectral shape of the non-thermal 
SZE as a function of the electrons spectral parameters. 
In this case, we use an arbitrary normalization in which 
the amplitude of the SZE at its minimum is set equal to 
— 1, and we vary the parameters of the electrons spectrum. 
Figl6]shows the SZE spectrum for a single power-law spec- 
trum 



fe{p) ocp 



P>Pi 



(30) 



(here p = f3"f is the normalized momentum), with si — 3.0, 
as a function of the minimum momentum pi . It is clear 
that the frequency position of the minimum of the SZE is 
different in the case pi = 1, while for pi > 5 the position 
of the minimum does not change for different values of pi . 
The shape of the spectrum at higher frequencies {ly > 300 
GHz) is, instead, more sensitive to the value of pi. 
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Fig. 6. The non-thermal SZE spectrum, normalized to its min- 
imum amplitude is shown as a function of the minimum mo- 
mentum pi of the electron spectrum. We assume here a single 
power-law spectrum with si — 3.0. 

We show in Figl7]the case of a double power-law spec- 
trum 



fe{p) 



P 



Pbreak (p/ Pbreak ) 



Pi <P< Pbreak 
"''^ p > Pbreak 



(31) 



with si = 0.1, Pbreak = 100 and S2 ~ 3.0. In this case the 
position of the minimum and the shape of the spectrum at 
frequencies ly < 400 GHz, depends weakly from the value 
of pi , while a stronger dependence from pi is present at 
higher frequencies (this dependence is anyway weaker than 
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for the single power-law case because of the lower power 
at the low-p end of the electron spectrum). 
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Fig. 7. The non-thermal SZE spectrum, normalized to its min- 
imum amplitude is shown as a function of the minimum mo- 
mentum pi of the electron spectrum. We assume here a double 
power-law spectrum with s\ = 0.1, s = 3.0 and pbreak = 100. 

Fig[5] shows the changes of the SZE spectrum in the 
case of a single power-law with pi = 1 as a function of 
the spectral index si: also in this case the shape of the 
SZE spectrum is more sensitive to the shape of the elec- 
tron spectrum at frequency i' > 400 GHz. Fig[9]show the 
variation of the SZE spectrum as a function of the spec- 
tral index S2 for a double power-law spectrum: in this case 
the dependence from the electron spectrum index is quite 
weak for the considered frequency range. 

After having shown the dependence of the SZE spec- 
trum on the electron spectrum parameters, we discuss now 
the procedures through which we can recover the physical 
parameters of the non-thermal electron spectrum from the 
SZE observations. We consider here two cases: i) the case 
of a pure non-thermal electron population with a power- 
law spectrum (a situation applicable to radio galaxy lobes, 
see e.g. Colafrancesco 2008b); ii) the case of a combi- 
nation of thermal plus non-thermal electron population 
with a double power-law spectrum (a situation applicable 
to galaxy clusters with non-thermal phenomena, see e.g. 
Colafrancesco et al. 2003). 

5.1. The case of a single power-law spectrum 

We describe here the procedure through which it is pos- 
sible to derive the parameters of an electron population 
with single power-law spectrum, assuming that this is the 
dominant electron population for the SZE production. 
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Fig. 8. The non-thermal SZE spectrum, normalized to its min- 
imum amplitude is shown as a function of the spectral index s\ 
of the electron spectrum. We assume here a single power-law 
spectrum with pi = 1. 
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Fig. 9. The non-thermal SZE spectrum, normalized to its min- 
imum amplitude is shown as a function of the spectral index S2 
of the electron spectrum. We assume here a double power-law 
spectrum with pi = 1, si = 0.1 and Pbreak = 100. 

Following the same procedure previously outlined for 
the thermal SZE in Sect. 4, we have sampled the non- 
thermal SZE produced at different projected radii and at 
six frequencies in the range 300 - 400 GHz where, as previ- 
ously noticed, the shape of the SZE is maximally sensitive 
to the physical parameters of the electron spectrum. The 
emission region has been assumed, for simplicity, to have 
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spherical symmetry with a radius of 50 kpc with an elec- 
tron spectrum 



J^e{p, r) = fcoP • ge{r) P>Pl 



(32) 



Here we assume si — 3.0, pi = 3, /cq = 10 ^ cm and a 
radial distribution of the electron population given by 



21 -1- 



(33) 



with Tc = 10 kpc and = 0.5. 

In the fitting procedure we discuss here, we fix the 
spectral index of the electron spectrum si, because it can 
be derived from radio observations, and we fit the simu- 
lated observations with two free parameters: the optical 
depth T and the minimum momentum pi. We note that 
the assumption of the parameter si from radio data is not 
strictly necessary, because, in the presence of adequate 
coverage in frequency, it can be derived from the SZE 
data alone; this assumption is made for the sole purpose 
of reducing the time of calculation. 

The SZE intensity at different radii and the relative sam- 
pling are shown in Fig llOl The behaviour of the optical 
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Fig. 11. The radial profile of the optical depth as obtained 
from the fit to the SZE spectra b using data sampled in the 
range 300-400 GHz with uncertainty of 0.1% (see FigfTO)l. 
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Fig. 10. The non-thermal SZE spectrum evaluated at various 
projected radii and frequency sampling (see text for details) 
with data uncertainty of 0.1%. Curves are, from bottom to 
top, for projected radii of 0, 5, 10, 15, 20, 30, 40 kpc. 

depth at different projected radii, as derived by the fit- 
ting procedure, is shown in Fig. [TT] Our analysis recovers 
the values pi — 3.1 ± 1.0 at all the considered projected 
radii. In order to deproject the optical depth radial distri- 
bution and derive the radial profile of the electron density 
it is sufficient to fit the optical depth profile with a density 
profile as in eq. ([55)1 and integrate along the line of sight. 
The deprojection procedure yields ko = (9.0 ±3.5) x 10^^ 



cm-3, rc = 7.2 ± 6.0 kpc and = 0.43 ± 0.21; these val- 
ues are consistent, within the errors, with the input values 
that we used in the simulations. 

We stress that the previous analysis of the SZE spec- 
tral measurements allows to derive directly the density 
profile and the spectral shape of the electron distribution 
in the radio lobes. Therefore, the combination of SZE spec- 
tral measurements and measurement of the low-frequency 
radio emission allows to break the degeneracy existing in 
synchrotron emission between the electron density and the 
magnetic field and eventually yield a measure of the mag- 
netic field in radio lobes (see also Colafrancesco 2008b). 

5.2. The case of thermal plus non-thermal populations 

Here we consider a combination of a thermal electron dis- 
tribution that fits the A2199 data and an additional non- 
thermal electron spectrum with a double power-law spec- 
trum, as in eq. ([3T|) . with parameters si = 0.1, S2 = 3.0, 

Pl = 1 and Pbreak = 100. 

A non-thermal electron spectrum similar to a double 
power-law is expected in galaxy clusters because the non- 
thermal electrons of relatively low energy loose energy 
through Coulomb interactions with the thermal IC gas, 
and these energy losses produce a flattening of the elec- 
tron spectrum at low ( ^ 150 MeV) energies. 
To study this case, we do not follow the same strategy 
described in the previous Sections: specifically, we do not 
perform a fit to the sampled points because this proce- 
dure should involve too many free parameters (i.e., those 
of both the thermal and the non-thermal electron popula- 
tions) . Instead, we focus on the spectral regions where the 
non-thermal effect is expected to be more relevant, and 
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therefore it can be measured. For this purpose, we com- 
pare the non-thermal SZE spectrum with the thermal SZE 
spectrum (calculated both at first order in r and with sec- 
ond order corrections, see Sect. 2 and Colafrancesco et al. 
2003 for more details), and with the spectrum of a possible 
SZ kinematic effect. The spectral regions where the non- 
thermal SZE could be better estimated are the frequency 
range around the minimum of the SZE (at ;y ~ 150 GHz), 
because all the amplitudes are basically proportional to 
the optical depth of the relative electronic population, and 
the frequency range around the crossover frequency of the 
thermal SZE (at v - 220 GHz), because the thermal SZE 
is by definition null at such frequency leaving hence visible 
the (negative in sign) non-thermal SZE. 

We note that an upper limit to the non-thermal elec- 
tron density can be derived, after fixing their spectral 
shape, by requiring that the non-thermal pressure does not 
exceeds the thermal one (see discussion in Colafrancesco et 
al. 2003). FigfT2l shows the ratio between the non-thermal 
SZE and the thermal SZE in the case Pnt ~ Pth at the 
cluster center: this case can be considered as an upper 
limit to the contribution of the non-thermal component 
to the total SZE. Fig. [13] shows the spectrum of the ther- 
mal and non-thermal SZE in the frequency ranges 100- 
200 GHz and 300-400 GHz. The ratio of the non-thermal 
to the thermal SZE is large only in the region around 
V ~ 220 GHz where the thermal SZE has its crossover 
and where, consequently, the total SZE is dominated by 
the non-thermal component. For v < 200 GHz the non- 
thermal SZE is ^ 1% of the thermal one, while for i' > 300 
GHz the ratio becomes very low, of order of ^ 0.1 — 0.5%, 
or even less. Two conclusions can be derived from these 
results: 1) In the frequency range 300-400 GHz, that is op- 
timal to derive precise information on the thermal SZE, 
the non-thermal component does not provide a contribu- 
tion that jeopardizes the precision of the measurements. 
2) At lower frequencies, around 200 GHz, the contribution 
of the non-thermal component is more relevant; however, 
to determine such a component, it is necessary to separate 
it from the thermal one. 

A possible observational strategy to separate the non- 
thermal effect from the thermal one could consist in mea- 
suring the thermal SZE in the frequency range 300-400 
GHz and deriving the IC gas physical parameters, as seen 
in Sect. 4. Then, once the thermal SZE parameters are 
fixed, it will be possible to measure the SZE in the fre- 
quency range 100-200 GHz and derive the parameters of 
the non-thermal population that, in this frequency range, 
yields a contribution to the total SZE of order of 1% 
(see figures [T2] and fT3|l : such SZE amplitude is measurable 
if measurements have an uncertainty of order of ~ 0.1%. 

For a precise determination of the cluster parameters, 
we must take into account the contribution of the thermal 
SZE evaluated at second order in r; this can be calcu- 
lated in details once the main parameters of the IC gas 
are known (see Colafrancesco et al. 2003 for details). The 
modifications induced by the second order effect, in com- 
parison to the non-thermal SZE and to the first-order ther- 
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Fig. 12. The frequency dependence of the ratio of the non- 
thermal SZE (we assume a double power-law spectrum with 
si = 0.1, S2 — 3.0, pi — 1 and ptreak ~ 100) to the thermal 
SZE at the center of A2199. We assume a pressure balance 
P„t ~ Pth at the cluster center. 

mal SZE, are shown in Fig. 1141 For the frequency range 
300-400 GHz, that is the optimal frequency window to ex- 
tract the thermal SZE parameters, the second-order cor- 
rection to the thermal SZE is less than 30% of the non- 
thermal SZE and it is less than 0.1% of the first-order 
thermal SZE. Hence, we can conclude that the total SZE 
can be adequately approximated, in this frequency range, 
by the first-order thermal SZE only. Moreover, in the fre- 
quency range 100-200 GHz, where the non-thermal effect 
could be measured, the second order correction to thermal 
effect is ~ 5% of the non-thermal effect, so that the second 
order correction can be neglected also for this purpose. 

We consider now also the possible contribution of a 
kinematic SZE whose amplitude is maximal at frequency 
around the crossover frequency of the thermal SZE, i.e. 
dX V ^ 220 GHz. Fig. [15] shows the ratio of the kinemat- 
ical to the non-thermal SZE (upper panel) and the ratio 
of the kinematical to the thermal SZE (lower panel) for a 
receding cluster with peculiar velocity of 1000 km/s. The 
result is that in the frequency range 300-400 GHz also 
the kinematical contribution is a small fraction 20%) 
of the non-thermal SZE, and it is less than 0.1% of the 
thermal one. Therefore, also the kinematical SZE can be 
neglected in this frequency range. 

Fig.[Tn]also shows that in the range 100-200 GHz the kine- 
matical SZE is less than 10% of the non-thermal SZE. 
However, if the non-thermal electron population have a 
pressure Pnt ^ 0-^Pth, then the non-thermal and kinemat- 
ical SZE would become comparable (under the assumption 
of having Vp = 1000 km/s); in such a case, however, these 
signals would be indistinguishable from the thermal one 
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Fig. 13. The thermal (solid curve) and non-thermal (dashed 
curve) SZE evaluated at the center of A2199. The non-thermal 
electron population has a double power-law spectrum with pa- 
rameters si = 0.1, S2 ~ 3.0, pi — 1 and Pbreak = 100 and 
a pressure P„t ~ Pth at the cluster center. The spectra are 
shown in the frequency range 100-200 GHz (upper panel) and 
300-400 GHz (lower panel). 

(see for reference. Figs. [T^] and [T3] upper panel, in which 
we show the case P„t ~ Pth)- 

We can summarize the results obtained: 

- in the frequency range 300-400 GHz, the first order ther- 
mal SZE is a good approximation of the total SZE (all the 
corrections are < 0.1%). 

- in the frequency range 100-200 GHz, if the non-thermal 
contribution is strong, it dominates over the kinematic 
one. If, instead, the two SZE contributions are compa- 
rable, then they cannot be easily disentangled from the 




V (GHz) 

0.0101 ' ' ' 'II 

0.008 - 

<f 0.006 - 




V (GHz) 

Fig. 14. Upper panel: the frequency dependence of the abso- 
lute value of the ratio between the second order correction to 
the thermal SZE at the center of A2199 and the non-thermal 
SZE of an electron population with double power-law spectrum 
with parameters si = 0.1, S2 = 3.0, pi = 1, Pbreak = 100 and 
with pressure Put ~ Pth- Lower panel: the frequency depen- 
dence of the ratio between the second order correction and the 
first order description of the thermal SZE. 



thermal one. Hence, if deviations from expected thermal 
SZE are measured in this frequency range, they should be 
of non-thermal origin. 

- only at ~ 220 GHz (i.e. the crossover frequency of ther- 
mal SZE) the kinematic SZE, the non-thermal one and 
the second order thermal correction are important with 
respect to the first order thermal effect. To disentangle 
between the three different contributions it is important 
to put constraints on the thermal and the non-thermal pa- 
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Fig. 15. Upper panel: the frequency dependence of the abso- 
lute value of the ratio between the kinematic SZE for a value 
of the peculiar velocity of Vp — 1000 km/s and the non-thermal 
SZE for an electron population with the same spectrum as in 
Fig. 1141 Lower panel: the frequency dependence of the absolute 
value of the ratio between the kinematic SZE and the thermal 
SZE evaluated at first order. 

rameters by performing precise observations in frequency 
ranges 300-400 and 100-200 GHz, respectively. 

6. Applications to specific clusters: Perseus and 
Ophiuchus 

We present in this section the results of our studies for 
two nearby and extended clusters: Perseus, which has a 
cool core, and Ophiuchus, which is approximately isother- 
mal. As we discussed in the previous sections, the optimal 
strategy to extract the physical parameters of the cluster 



is to consider the frequency range 300-400 GHz; in this 
range, we can safely consider only the first order thermal 
SZE in our analysis, because higher-order corrections are 
negligible. 

6.1. The Perseus cluster 

The Perseus cluster is at a redshift z = 0.0179 or at a 
distance of 77.7 Mpc, at which 1 arcmin corresponds to 
21.8 kpc. Therefore, a study of the spatially resolved SZE 
allows to have a spatial sampling of the temperature pro- 
file better than in the case of A2199 previously discussed. 
In the following, we repeat the previous analysis for this 
cluster. 

Churazov et al. (2003) derived an analytical approxi- 
mation for the radial profile of the IC gas density and of 
the gas temperature that fits the cluster data for r > 10 
kpc, while at shorter distances the presence of the central 
dominating galaxy NGC1275 makes the IC gas radial pro- 
file more uncertain. Here, for sake of illustration, we use 
the analytical fitting formulae given by Churazov et al. 
(2003) and we estrapolate them down to the cluster cen- 
ter, e.g. ignoring, in this way, the presence of the central 
galaxy. 

We have hence simulated the thermal SZE to first or- 
der expected from Perseus at seven projected radii from 
the cluster center: 0, 20, 50, 75, 100, 150 and 200 kpc 
(corresponding to angular sizes of 0, 0.9, 2.3, 3.4, 4.7, 6.9 
and 9.2 arcmin). We have then sampled each one of the 
SZE spectra at these radii for eight different frequencies: 
300, 320, 340, 360, 380, 400, 420 and 440 GHz (see Fig. 
I16|) . We have extended the frequency coverage in order to 
better sample the spectral region of the SZE where rela- 
tivistic effects are relevant. For each projected radius we 
have finally fitted, using the gas temperature and opti- 
cal depth as free parameters, the SZE simulated signal by 
using the expression of the first-order thermal SZE and 
assuming an error of 1% in the experimental data. The 
results of the fit for the projected temperature are shown 
in Fig. [T71 The deprojection procedure of the best-fit ra- 
dial temperature profile is done by using Eq. ([28]) setting 
the values kBT^xt — 6.8 keV (the best-fit projected tem- 
perature at 200 kpc) , and by using as prior constraints the 
information on electrons radial profile provided by X-rays 
observations (see Sect. 4); for an experimental uncertainty 
of 1%, such procedure yields values ksTint = 3.6±1.5 keV, 
rc = 104±45 kpc, = 4.0±2.4, with = 0.1 for 4 d.o.f. 
Fig. [m shows the comparison between the deprojected 
temperature profile found with our method and the best- 
fit deprojected temperature profile derived by Churazov 
et al. (2003). The two curves are quite similar except for 
the inner region where the IC gas is cooler and the sensi- 
tivity of the SZE spectrum to low temperatures (through 
the effect of relativistic corrections) is lower. We have also 
verified that assuming an uncertainty of 0.1% in the ex- 
perimental data, the best-fit values do not change, while 
the relative errors are considerably reduced (see Fig. [19]) . 
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Fig. 16. The thermal SZE spectra of the Perseus cluster eval- 
uated at different projected radii of 0, 20, 50, 75, 100, 150 
e 200 kpc (from top to bottom) are shown in the frequency 
range 200-500 GHz. For each curve the SZE sampling in the 
frequency range 300-450 GHz has an uncertainty of 1%. 



Fig. 18. The deprojected temperature radial profile in Perseus 
derived from the fit to the X-ray data by Chrurazov et al. 
(2003) (solid curve) is compared to the fit to the SZE simulated 
data (dashed curve) shown in Fig.[l7l For the SZE temperature 
profile we use the radial profile given by eq. (|28l) . with ksTcxt ~ 
6.8 keV, fcsTint = 3.6 keV, rc = 104 kpc and /i = 4.0. 
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Fig. 17. The radial profile of the projected temperature, as 
derived from the fit to the SZE spectra in Perseus in the fre- 
quency range 300-440 GHz with uncertainties of 1% (see Fig. 
I16p . is shown as a function of the projected radius (in arcmin). 



6.2. The Ophiuchus cluster 

The Ophiuchus cluster is, at first approximation, an 
isothermal cluster with kT ~ 9.9 keV (Watanabe et al. 
2001). Moreover, being Ophiuchus a quite hot cluster, rel- 
ativistic effects in the SZE spectrum are expected to be 



Fig. 19. As Fig. [17] but with uncertainties assumed at the 
level of 0.1% on the sampled SZE data points (see Fig. I16[l. 

more prominent and this allows to obtain a better esti- 
mate of both the cluster temperature and density. 

The SZE signal produced in Ophiuchus, calculated as 
described in Sect. 2, has been sampled at six frequencies 
between 300 and 400 GHz at seven different radii: 0, 35, 70, 
120, 200, 300 and 500 kpc (see FigEOl), assuming experi- 
mental uncertainties of 1%. The minimum radius at which 
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we sampled the radial profiles has been chosen taking into 
account that at the Ophiuchus distance (z = 0.028), 1 
arcmin corresponds to ^ 34 kpc. 









Fig. 20. The thermal SZE spectrum for the Ophiuchus cluster 
is shown at various projected radii of 0, 35, 70, 120, 200, 300 
and 500 kpc (from top to bottom). The assumed uncertainties 
on the SZE data points is of 1%. 



FigHH shows the radial profiles of the projected tem- 
perature and optical depth obtained from the fit to SZE 
data. These results coherently indicate that the cluster 
temperature is constant with radius, within the errors, 
while the radial density profile can be derived from the 
optical depth radial profile by using the following model 
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Therefore, we note that in a cluster with constant tem- 
perature profile, the deprojection procedure can be made 
by using only optical depth results without assuming any 
prior constraint taken from X-ray data. Hence, the ther- 
mal gas parameters derived from SZE and X-ray data can 
be directly compared. 

The values derived from X-ray measurements are 
rieo = 1.77 X 10-2 cm-3, Tc = 108 kpc and qth = 0.96 
(Johnston et al. 1981; Watanabe et al. 2001). The values 
we derive from the fit to the simulated SZE signals are 
rieo = (1.7 ± 0.2) X 10-2 cm-3, = 118 ± 2 kpc and 
qth = 1-00 ± 0.07. The comparison between the best-fit 
profile to X-ray data and the one to the SZE data is shown 
in Fig[2H We can see that the density profile derived from 
SZE measurements is very similar to that obtained from 
X-ray observations. 



Fig. 21. The radial profile of the projected temperature (up- 
per panel) and of the IC gas density (lower panel) as derived 
from the fit to the SZE spectra for Ophiuchus is shown as a 
function of the projected radius (in arcmin). We use SZE data 
points in the frequency range 300-400 GHz with uncertainties 
of 1% (see Fig. [Mil. 

7. Discussion and conclusions 

We derived in this paper a general formalism for describing 
the thermal SZE in galaxy clusters with a non-uniform 
temperature profile that can be applied to both cool-core 
clusters and non-cool core cluster with an isothermal or 
non-isothermal temperature structure. 

In order to extract the temperature profile information 
from spatially resolved, spectroscopic SZE observations, 
we derive an inversion technique through which the elec- 
tron distribution function can be extracted using a wide 
frequency range coverage of the SZE signal. We present 
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Fig. 22. The IC gas density profile of Ophiuchus as derived 
from X-ray observations (solid curve) is compared with the one 
derived from the optical depth radial profile shown in Fig |21l 
(dashed line) obtained from SZE simulated observations (see 
text for details). 



extensive results of the fitting procedure used to extract 
the cluster temperature from a set of simulated spatially 
resolved spectroscopic SZE observations in different bands 
of the spectrum, from 100 to 450 GHz, and we discuss the 
optimal experimental and theoretical strategy. 

The results of our analysis for three different clus- 
ter prototypes (A2199 with a low-temperature cool core, 
Perseus with a relatively high-temperature cool core, 
Ophiuchus with an isothermal temperature distribution) 
provide the required precision of the SZE observations and 
the optimal frequency bands for a determination of the 
cluster temperature similar or better than that obtain- 
able from X-ray observations. The precision level of SZE- 
derived temperature is also discussed for the outer regions 
of clusters, an issue that is particularly relevant for accu- 
rate mass determination of galaxy clusters. 
We also study the possibility to extract from our method 
the parameters characterizing the non-thermal SZE spec- 
trum of the relativistic plasma contained in various astro- 
physical environments: the lobes of radio galaxies as well 
as the relativistic electrons co-spatially distributed with 
the thermal plasma in galaxy clusters with non-thermal 
phenomena (e.g. radio halos and/or hard X-ray excesses). 

Based on our results, we conclude that the method 
we present here is the only method, so far, for SZE ob- 
servations that is able to extract the crucial parameters 
of the cluster atmospheres (i.e. their temperature, density 
and additional non-thermal components) using only a sin- 
gle observational technique, i.e. spatially resolved spectro- 
scopic SZE observations. 



Other methods to extract the cluster temperature 
profile by using SZE observations have been presented 
(Holder & Loeb 2004), but they make use of a combination 
of the thermal SZE (which measures the electron pressure 
distribution) and of radio observations of the highly po- 
larized scattered radiation coming from active galaxies as- 
sociated to the cluster (which depends on the cluster elec- 
tron density distribution) . Holder & Loeb pointed out that 
current instruments should allow us to reach accuracy of 
the mass- weighted cluster temperature profiles of order of 
~ 1 keV, but under the assumption that the central radio 
source is steady over several million years. However, vari- 
able or beamed sources (like radio galaxies or jets/lobes of 
AGNs associated to the cluster) will leave observable sig- 
natures in the scattered emission. Therefore, this method 
could allow, in principle (e.g. with fiK sensitivity obser- 
vations to polarized emission) , to measure the age of the 
central source by finding an edge to the polarized emis- 
sion. 

We noticed in this context, that combining this method 
with the independent temperature measurements from 
SZE spectral studies (that we discuss here) would allow 
to measure the time evolution of the central radio source. 

In conclusion, our study shows that the next 
generation SZE experiments with spectroscopic 
capabilities, like those using FTS spectrome- 
ters with imaging capabilities (e.g. SAGACE, 
\http://oheron.romal. infn.it/sagacemission/ , and/or 
MILLIMETRON), can provide precise temperature 
distribution over large radial distances for galaxy clusters 
even out to substantial redshift. This will allow to use 
SZE observations of clusters to better understand the 
physics of cluster atmospheres, to better reconstruct their 
total mass content and, eventually, to use clusters of 
galaxies as reliable astrophysical and cosmological probes. 
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Appendix A: Comparing our results to the 
covariant formalism derivation of the SZE. 

We compare here the SZE calculated according to the for- 
malism of Wright (1979; W79) with the one calculated 
according to the derivation of Boehm & Lavalle (2009; 
BL09). For this last case, we refer to their eqs. (18)-(22), 
derived in the Thomson limit, that are therefore directly 
comparable with those derived by W79. 

According to BL09, the correct expression of the SZE 
in the relativistic covariant formulation, is: 



(A.l) 



The averages over the angles of the two contributions are 
given by the following expressions 



/^'(^^fc) =2/CTnr/°(£;fc), 

where /C ^ 1 in the Thomson limit, and 



(A.2) 



2t,,, j dpfe{Ep) 



(Tt 



(A.3) 



Here, the function fe{Ep) is normalized as to give 
Jd^p/{27rffe{Ep) = 1, where t = (1 ~ /3/i)/(l - /3/i'), 
and the function is given by: 



J^(/3,/i,A*') 



f3^m^ 3aT (1 - Pn') 
(27r)3 16 (l-/3Ai)2 

2 - 2K{1 - ^1^') + if 2 [(1 - 



-^(i-^^)(i-m'^) 



(A.4) 



where K = [7^(1 - ^^)(1 - |3^l')]-^. 

For the quantitative comparison of the two approaches, 
we have used the case of an electron population with a 
double power-law spectrum: 



fe{p) 



— Si 
Pbreak 



{p/ Pbreak 



Pl<P< Pbreak 
"■'^ p > Pbreak 



where p = f3j, that is normalized as to give Jp°° fe{p)dp = 
1. We use parameters si = 0.1, S2 = 3.0, pi = 1 e 

Pbreak = 100. 
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We notice that, with the use of the a-dimensional mo- 
mentum p and the previous normahzation, the expres- 
sion for the function in eq. (|A.4p shghtly 
changes: for consistency it must be multiplied by the fac- 
tor (27r)^/(/3^m^7^), as derived by using eqs. (16), (17) 
and (53) in BL09. 

Fig. lA.ll shows the SZE spectrum plotted over a wide 
frequency range, x — {hv) / {ksTQ) from x = \Q [v ^ 570 
GHz) up to cc = 10^5 {hv - 2.35 x 10^ GeV) calculated 
up to first order for t — 1 according to W79 (continu- 
ous curve) and that calculated according to BL09 (dashed 
line). As is clearly recognized, the two spectra are indis- 
tinguishable over the whole frequency range. 




Fig. A.l. The SZE spectrum calcolated following the method 
of W79 (solid line) and following the method of BL09 (dashed 
line). The horizontal axis shows the a-dimensional frequency 
X = (hu) / [kBTo) while the vertical axis shows the SZE bright- 
ness change Ai(x) in units of 2(fcsTo)^/(/ic)^- No difference 
between the two approaches appears. 

Fig. IA.2I shows a blow-up of the SZE spectra in the 
microwave region (x = 0—20). Also in this region we must 
conclude that the results obtained with the two formalisms 
are indistinguishable. 

Fig. I A. 3 1 shows the comparison between the ICS emis- 
sion (against the CMB photons) along the line of sight 
calculated in the approach of BL09 and in the standard ap- 
proach (see e.g. Colafrancesco et al. 2005, see also Longair 
1993) using also the Klein-Nishina cross-section, and as- 
suming an electron spectrum as in eq. (jA.5|) , with a density 
and spatial distribution such that t — 1. The two dif- 
ferent calculations provide quite similar results, apart for 
the frequency range x > 10^^, at which the Klein-Nishina 
cross-section has a non-negligible effect (it decreases, in 
fact, the Compton scattering efficiency at high energies, 
see e.g. Fargion & Salis 1998). 



0.5 




-1.5 1 

5 10 15 20 



X 



Fig. A. 2. Same as Fig. lA.ll but with a zoom in the microwave 
frequency range at which the SZE telescope operate. Also in 
this low-frequency range, there is no difference between the two 
approaches. 




Fig. A. 3. The comparison of the ICS emission evaluated fol- 
lowing the standard computational approach with the inclu- 
sion of the Klein-Nishina cross-section (solid line) and the one 
evaluated following the method of BL09 in the Thomson ap- 
proximation (dashed line). No difference appears in the two 
approaches up to a frequency x ~ lO'^^. 

Based on the previous results, we can conclude that: 
i) for the SZE spectrum in the microwave-mm region the 
formalism of W79 can be used also in the case of relativis- 
tic electrons since it provides the same results of the co- 
variant formalisms of BL09. In addition, the W79 formal- 
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ism, that uses a more simple expression requires a much 
shorter computing time w.r.t. the BL09 formalism, 
ii) for the calculation of the ICS emission at high frequen- 
cies, the standard formalism is completely adequate since 
also in this case there is no difference w.r.t. the covariant 
formalism results. In addition, the standard formalism re- 
quires much shorter computing time and allows also the 
use of the Klein- Nishina cross-section effect without intro- 
ducing additional complications in the numerical compu- 
tation, as occurs - on the contrary - in eq.(46) of BL09. 

Recently Nozawa & Kohyama (2009) analyzed the co- 
variant formalism of the Sunyaev-Zeldovich effect for the 
thermal and nonthermal distributions and derived the 
frequency redistribution function identical to the W79 
method assuming the smallness of the photon energy (in 
the Thomson limit). They also derive the redistribution 
function in the covariant formalism in the Thomson limit. 
These authors have shown that two redistribution func- 
tions are mathematically equivalent in the Thomson limit 
which is fully valid for the cosmic microwave background 
photon energies. 

To summarize, the W79 and the covariant formalisms 
to calculate the SZE arc fully equivalent, contrary to pre- 
vious erroneous claims (see BL09), and the large advan- 
tage given by the standard calculations based on the W79 
approach are due to the much shorter computing times. 



